Merging of Dirac points and Floquet topological transitions in AC driven graphene 
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In the last years, the search for new topological phases in condensed matter sys- 
tems has become a priority. In particular, their control by external means is still an 
important challenge from both fundamental and technological point of view. On the 
other hand, graphene, because of its striking electronic band structure, is an excellent 
platform to explore such phases. In particular, the manipulation of its band structure 
is an important issue. We demonstrate that driving graphene w^ith an in-plane ac 
electric field allow^s for tuning the band structure and then induce out-of-equilibrium 
semi-metallic and insulating phases with distinct topological properties. In particular, 
suitable drivings offer the possibility to move and merge the Dirac points of graphene 
— a topological transition which has not been observed yet in electronic systems — but 
also to localize the electrons in specific directions, or even to increase the number of 
Dirac points. 

PACS numbers: 



Introduction - Since the recent discovery of topo- 
logical insulators [1, 2j, the search for topological transi- 
tions in condensed matter has become a priority task. In 
particular, the prediction and the observation of a new 
topological order, namely the Z2 topological order [3-5J, 
has stimulated a great interest in the scientific commu- 
nity. Graphene, a two-dimensional carbon-based crys- 
tal, is a wonderful platform for investigating topologi- 
cal transitions. It is a semi- metal, whose band structure 
of the Pz orbitals consists in two bands that touch lin- 
early at the Fermi level at two inequivalent points of the 
Brillouin zone [6j. This peculiar structure gives rise to 
massless Dirac-like low energy excitations. The quan- 
tum spin Hall effect, reflecting the Z2 topological order, 
was actually first proposed in graphene by Kane and 
Mele [7J , but the strength of the intrinsic spin-orbit cou- 
pling is unfortunately too weak to induce a reasonably 
large inverted gap at the Dirac points, where the bands 
touch, and lead the system to a topological insulating 
phase. Recently, the idea of inducing the band inversion 
in HgTe/CdTe heterostructures by means of an exter- 
nal potential, periodic in time, rather than by changing 
the width of the quantum well, was proposed by Lin- 
der et al. [8j. Such an out-of-equilibrium topological 
phase has been recently observed in temporal modulated 
photonic crystals [9j. This approach renews considerably 
the possibilities of inducing new topological phases. An 
interesting application is of course graphene, where topo- 
logically non-trivial insulating phases may be driven by 
the interaction between the electron charge and the AC 
electric field [10, llj. Remarkably in graphene, a differ- 
ent topological transition between a Dirac semi-metallic 
(SM) phase and an insulating (I) phase, can be achieved 
by merging the two Dirac points p^ [13j . The resulting 
I phase is not a Z2 topological phase, but may neverthe- 
less hosts zero-energy edge states whose topological ori- 



gin is well understood in terms of a one-dimensional bulk 
winding number, namely the Zak phase [14HT6] . Such 
a topological merging transition is reached in a honey- 
comb lattice whenever a critical homogeneous anisotropy 
is introduced. The search for this transition has stimu- 
lated experimental efforts beyond the solid state commu- 
nity and the merging and alignement of Dirac points has 
been theoretically analyzed in shaken honeycomb opti- 
cal lattices for the high frequency regime[17j. Also the 
merging of the Dirac points as well as the emergence of 
zero-energy edge states have been recently respectively 
observed in anisotropic traps of cold atoms [T8[ [T9] and in 
microwave tight-binding analogue experiments of an hon- 
eycomb lattice [20IEI]. However, the merging transition 
has not been observed yet in electronic systems. In par- 
ticular its achievement in graphene itself [22j, by means 
of mechanical manipulations like stretching, is unfortu- 
nately out of reach, because the graphene sheet would be 
destroyed far before the expected transition [23j. 




Figure 1: Graphene sheet irradiated by an AC electric field 
with arbitrary polarization. The inset shows the merging of 
the Dirac points induced by the external field. 



In the present work, we show that an AC electric 
field can induce the topological merging transition 
in graphene. More precisely, we demonstrate that in 
the high frequency regime, the AC field acts similarly 
to a mechanical strain, and therefore allows for the 
manipulation of the Dirac points, their annihilation 
and their creation, in a controllable way. Furthermore, 
the AC field is able to induce the localization of the 
electrons in specific directions, which marks another 
topological transition between two I or two SM phases. 
Our analysis for the high frequency regime with lin- 
ear polarization agrees with the results obtained in 
shaken optical lattices [17J. In addition, we also consider 
different field polarizations. As we will see, it plays 
a crucial role for lower frequencies by breaking time 
reversal symmetry (TRS). The analytical expressions 
obtained here provide a further insight out of the high 
frequency regime. At lower frequencies, we find that the 
coupling between the Floquet bands, that characterize 
the properties of the driven system, leads to Floquet 
bands inversions. This gives rise to out-of-equilibrium 
phases with no analog in the static system in which 
e.g., multiple PDF appear. Importantly, our approach 
does not rely neither on low energy models nor on 
approximations valid just close to resonance (as the 
rotating wave approximation) [8l \TT\ [24j . In particular, it 
allows for keeping track of the relative position between 
the Dirac points due to the interplay between lattice 
and field symmetries. In addition, the results are valid 
for arbitrary field amplitudes and phase polarizations. 

Floquet theory of graphene - The Hamiltonian of 

a T-periodic driven system fulfills H{r-\-T) = H (r), 
and its eigenvectors can then be written in the Flo- 
quet form: \^p {r)) = e~'^^^\^ (r)) {h = 1), being e the 
quasi-energy, and |<^ (r)) = \^ {r -\-T)) the Floquet state 
|25j. This ansatz, consequence of the time translation 
invariance, maps the time-dependent Schrdinger equa- 
tion to the eigenvalue equation H (t) |$ (r)) = e|$ (r)), 
in which the quasi-energies and the Floquet states are 
the eigenvalues and eigenvectors of the Floquet operator 
H{r) = H (r) — idr respectively. In order to deal with 
the time-dependence of the Floquet operator H (r), it is 
useful to introduce the Sambe space, which considers the 
Floquet states as belonging to a composed Hilbert space, 
where the basis states are time-independent [26j. In this 
space, the quasi-energies are given by (see supplementary 
materials): 






/puj 



(1) 



where Hp^-p 



f^ ^e'^^^P-P^H{T) with CO = 27r/T, and 
where ^pio is the p^'"^ — th G Z Fourier component of 
the Floquet state. Eq. [T] shows that the quasi-energies 
structure consists in copies of e G [~f 7 f]? oi" Floquet 



side-bands, shifted by pu, which reflects the periodicity 
of the driving. The coupling between the Floquet side- 
bands is encoded into the Fourier components Hp>-p of 
the Hamiltonian. 

We now consider the concrete case of a graphene sheet 
embedded in an in-plane periodic time-dependent electric 
field f (r), as displayed in figure [l] The crystallographic 
structure of graphene consists in a triangular Bravais lat- 
tice with two atoms per unit cell separated by a distance 
a. In the basis of the two sub-lattices, the time-dependent 
tight-binding Bloch Hamiltonian has the general form: 



/J(T,k) 



p(r,k) 
p*(r,k) 



(2) 



The wave function of an electron coupled to an elec- 
tric field £ {t) = —drA{r) carries a phase f d£ • A(r) 
along the path d£. Thus, the hoppings t to the three 
nearest neighbors on the undriven honeycomb lattice be- 
come t -^ tj{^) = te'^^^'-^^^'^ where the vectors joining 
the nearest neighbors dj are given by di = | (—1, a/3), 
d2 = f (-1,-^3), da = a (1,0). The /)(r, k) function 
entering the Bloch Hamiltonian therefore reads: 



p(r,k) =ti{r)e 



;kai 



t2{T)e 



ika2 



hir) 



(3) 



with ai and a2 the basis vectors of the Bravais lat- 
tice. We consider electric fields with arbitrary ampli- 
tude and polarization, and write the vector potential as 
A(r) = {Ax sm{ujr) ^ Ay sm{ujr -\- cp)) . The analytical cal- 
culation of the Fourier components Hq^pf-pCk) of the 
Bloch Hamiltonian leads to (see supplementary materi- 
als): 



Hg{k) = 



p,(k) 
P-,(k) 



(4) 



5]] • t^ e*^'^^ , where we have defined the 



with pq{k) ^^ .^^q 

time independent Floquet hoppings as: 

7 q — iJ — qy^jjC 



(5) 



Jq (x) is the q^^ order Bessel function of the first kind, and 
the dimensionless functions Aj and ^j which encode all 
the information about the amplitude and the polarization 
of the field are defined as: 



A2,l 
^2,1 

As 



Al + 3A2 ± 2V?>AxAy cos ((/:)), 



= ± arctan 



y/?>Ay sin ((/:?) 
Ax ± V^Ay COS ((/?) 
= . 



(6) 



As a consequence, spatial anisotropy can be tuned in the 
driven system by varying the polarization or the ampli- 
tude of the field. This is the key result of the present 
work. 



For the sake of clarity, we first focus on the high fre- 
quency regime uo ^ t. In that limit, the Floquet bands 
are uncoupled and the zeroth Fourier component Hg=o of 
the Floquet operator is the dominant one. Thus, Eq{T]is 
block-diagonal in the Fourier space, and simply consists 
in a collection of 2 x 2 time-independent Hamiltonians 
separated in energy by uj. The quasi-energy of a side- 
band Of, is then simply given by: 

ea{\<i]A^,Ay,ip) = ±\pq=o{^;Ax^Ay,(p)\ ^au. (7) 

This shows that the quasi-energy of each Floquet 
band of AC driven graphene is, in the high-frequency 
regime, identical to the energy bands of undriven 
graphene with renormalized hopping parameters 

tf^Q{A^,Ay,ip) =tJo{Aj). 

Manipulation of the Dirac points by the AC 
field: merging and localization - The band structure 
of undriven graphene consists in two bands that touch 
linearly at two inequivalent points K^'"^ located at the 
corner of the Brillouin zone, the Dirac points. A direct 
consequence of the anisotropy of the hopping parame- 
ters induced by the AC field is to change the location 
of these Dirac points. The two Dirac points, related by 
time-reversal symmetry - and thus located at (/c, — /c) in 
the Brillouin zone - will move as the anisotropy is mod- 
ified, and merge at one of the four time-reversal sym- 
metric points of the Brillouin zone M^, whenever a spe- 
cific relation between the hopping parameters is fulfilled 
[131221 [27]: 



Mo, tf + tf +tf = 



Ml, tf = tf+tf 



Mo 



t2 — ti -\- ^3 



M, 



^3 — ^1 + ^2 



(8) 



where the index has been dropped out for clar- 
ity. The merging transition corresponds to the cre- 
ation/annihilation of a pair of Dirac points (PDF). At 
the transition, the dispersion relation is quadratic in one 
direction, but remains linear in the other one, leading to 
the prediction of striking properties such as an unusual 
temperature dependence of the specific heat and mag- 
netic field dependence of the Landau levels [27J. Such 
transitions can now be achieved for specific values of am- 
plitude {Ax, Ay) and polarization (p of the electric field. 
AC driven graphene exhibits several distinct SM and 
I phases that we now describe. Figure |2] shows a typical 
example of a SM/I phase diagram obtained in the high 
frequency regime when varying the amplitude and the 
polarization of the field (see supplementary materials for 
other exemples). The creation/annihilation of a PDF at 
a Mi point is represented by a continuous merging line 
that separates a SM phase from an I phase. There are 
therefore four different merging lines, one for each point 
M^. It follows that distinct SM or I phases may emerge 
when different merging lines are crossed. This distinction 



can be made in terms of both the points M^ and topolog- 
ical (winding) numbers. Let us first precise this point for 
the SM phases: Each SM phase can actually be labelled 
by a point M^ where a PDF cannot be annihilated. For 
instance, a PDF cannot be annihilated at Mq (green line) 
from the SM phase SMq (which is the one connected to 
the undriven graphene phase) . The reason is that all the 
hopping parameters in the SMq phase have the same sign, 
so that the merging at Mq would require first a change in 
sign of any of them, which, according to Eq. [8J would ful- 
fill the merging conditions at other M^^o points. There- 
fore, one can distinguish four SM phases, denoted as SM^, 
one for each M^ point where the merging transition can- 
not occur. In addition to this first criteria, two winding 
numbers can be introduced to characterize the topology 
of these four SM phases. The first one is the quantized 
Berry phase evaluated around a Dirac point, giving rise 
to a topological charge ^ f dk • Vk^k = ±1 assigned to 
each Dirac point, where 0\^ = argp(k) is the polar angle 
parametrizing the Bloch sphere. The two Dirac points of 
graphene carry opposite topological charges that annihi- 
late at the merging transition [22]. Besides, the topolog- 
ical transition is accompanied with a change of a second 
winding number, the Berry phase evaluated across a re- 
duced one- dimensional Brillouin zone, namely the Zak 
phase Z = ^ /_g/2 <^k * Vk^^k, where G is a vector of 
the reciprocal lattice. Chiral symmetry guaranties inte- 
ger values of the Zak phase [15j which depend on the 
direction in the reciprocal space. This reflects the edge 
orientations dependence for the density of zero-energy 
edge modes [16j. Thus, the four SM phases differ by the 
set of values the Zak phase takes for all directions in the 
reciprocal space, and therefore by the range of existence 
in /c-space of edge zero-energy modes for a given edge 
orientation. 

Similarly, different I phases can be distinguished as well 
by the set of values the Zak phase takes in all possible di- 
rections, even though a PDF can, in principle, be created 
at any of the four points M^. We find one normal insu- 
lating (NI) phase for which the Zak phase is zero in ev- 
ery directions (meaning the absence of zero-energy edge 
states), and, for polarizations different than (p = 7r/2, 
two Zak insulating (ZI) phases for which the Zak phase 
is non zero in different directions. The ZI phases reflect, 
in 2D, the underlying non-trivial topology expected for a 
BDI symmetry class in ID [28]. 

Interestingly, the phase diagrams show multiple cross- 
ings between the merging lines, meaning that a PDF can 
be created and annihilated simultaneously at two different 
Mi points. Following Eq. [8J such crossings actually im- 
ply the vanishing of one of the three hopping parameters 
tj, what directly induces localization of the electrons in 
the direction dj. This striking feature corresponds to a 
transition between two I phases or two SM phases. At 
the transition, the gap closes along lines parallel to the 
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Figure 2: SM/I phase diagram ((/?, aAy) in the high frequency 
regime for fixed A^jAy — v^, with a the lattice spacing. The 
different phases are separated by four distinct merging hues 
corresponding to the four M^ points of the first Brillouin zone 
(shown on the right) where a PDF can be created/annihilated. 



dj direction that passes through the two distinct merg- 
ing points. Unlike a single merging transition, it follows 
that the dispersion relation is flat along the dj direction 
but remains linear in the other one, as shown in Fig|3] 
(c). At the transition, the system is then reduced to an 
array of uncoupled one-dimensional chains, hosting one- 
dimensional massless Dirac fermions. Remarkably, for 
phase polarization Lp = 7r/2, the merging at Mi and M2 
always coincide. This gives rise to critical localization 
lines that separate two SM phases as shown in Fig. |3] 
At the transition between two SM phases, the topologi- 
cal charges assigned to the Dirac points change sign (see 
FigJ3](b) and (d)). In that sense, the localization (or dou- 
ble merging) transition is also a topological transition. 
Finally, we report that the phase diagram for polarization 
(p = 7r/2 also exhibits crossings between the four merging 
lines simultaneously[3_3j. At such critical points, all the 
hopping parameters tj vanish and the quasi-energy of 
each side-band is perfectly flat which leads to full charge 
localization in all directions. 



Multi side-bands effects and topological transi- 
tions - Up to here, the high frequency regime allowed 
us to neglect the couplings between the Floquet side- 
bands. Starting from this limit, we now comment about 
the change of the phase diagrams when the frequency is 
decreased. 

The high frequency regime holds as long as the fre- 
quency is larger than the bandwidth of the undriven sys- 
tem, namely, uj > 6t. For lower frequencies {3t < u < 6t), 
the Floquet band a = overlaps the two nearest side- 
bands a = ±1, as schematically illustrated in Figj4j In 
that case, the coupling between the Floquet side-bands 
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Figure 3: (a) SM/I phase diagram in the high frequency 
regime for Lp — 11/2^ where only three phases appear, due 
to the overlap between the Mi and M2 merging lines, (c) 
This overlap gives rise to localization lines separating two SM 
phases. At the transition between these two SM phases, the 
topological charges assigned to each Dirac point change sign 
((b) and (d)). This topological change at the transition goes 
together with a change of the value of the Zak phase, which 
is 1 (0) in the SM3 (SMq) phase along the path represented 
by a full (dashed) arrow. 



becomes relevant, and the simple high frequency effective 
theory is not accurate anymore. Besides, unless the po- 
larization of the field is linear, these couplings will break 
time reversal symmetry and will therefore make the ex- 
istence of Dirac SM phases much more involved [34J . In- 
deed, for phase polarization (^ ^ 0, the SM phases are 
perturbed such that a small gap opens. However, simi- 
larly to the Boron Nitride model, the topological charges 
assigned to the non-trivial Berry phases in each valley 
remain, so that their annihilation/creation can still be 
achieved by varying the parameters of the field. The 
resulting phase diagrams reveal only small qualitative 
changes for uo > 3t, in comparison with those obtained 
in the high frequency regime. In what follows, we will 
focus on the Dirac SM phases. For that reason, we shall 
consider linearly polarized fields only. 

The overlap between the Floquet side-bands allows for 
new hoppings which can now include the absorption or 
the emission of photons. These new processes produce 
transitions mediated by photons from the valence to the 
conduction band. At a; = 3t, the side-bands a = -\-l and 
a = —1 touch at e = for vanishing amplitude of the field 
(FigJ4J center), thus giving rise to a new crossing at the 
Mo point. This new crossing corresponds to a band in- 
version, in which the symmetry of the conduction and va- 




Figure 4: Schematic figure for the quasi-energies overlap 
considering the side-bands a — — 1,0,+1, as a function of 
the field amplitude A for three different frequencies. The 
side-band a = is enhanced by a colored area {[—uj/2]uj/2]). 
For uj — At (left), the nearest neighboring side-bands weakly 
couple, slightly modifying the merging phase diagram. For 
cj = 3t (center) a band crossing between side-bands a — -\-l 
and a = — 1 happens at e = 0, leading to a band inversion 
that changes the I or SM state of the system. For uj — 2t 
(right) we show in blue color the inverted band region. 



lence band is exchanged. In this case, the high frequency 
approximation fails, and the side-band inversions must 
be taken into account. The quasi-energy band structure 
in the first Floquet-Brillouin zone (FBZ) after this side- 
band inversion at the Mq point has occured is displayed 
in Figj5](left). This phase reveals the emergence of a new 
PDF related with photon mediated transitions. More in- 
teresting: a further decrease of the frequency leads to a 
second crossing between the side-bands at the M3 point 
which allows one to create a second PDP, also associated 
with a photon mediated transition, as shown in Fig|5] 
(right). 




Figure 5: Left: Emergence of new PDF at the Mo point in 
the first FBZ when the frequency is decreased to cj = 2. It 
from a NI phase in the high frequency limit, with a Ay = 3.6, 
Ax = and cp = 0. Right: Appearance of a second PDF 
at the M3 point for uj — 1.5t. These two new PDP both 
originate from the touching between the side-bands a = +1 
and a = — 1, which is related to photon mediated transitions 
(see fig[4|. 



Experimental realization - The high frequency 
regime {uj ^ t) is the most appropriate one to control the 
merging of the PDP. It requires field frequencies in the 
near ultraviolet. At these frequencies, the effect of the 
magnetic field associated to the electromagnetic radia- 



tion can reasonably be neglected in micron size samples. 
This frequency range is achievable, as well as the field 
amplitudes required to cross a merging line. Those are 
given by the condition q£oa/3t > 1 ^ £^0 > 3.41VA , 
being Sq the electric field amplitude and a = 2.46Ap9l- 
I3T] . On the other hand, even though a heating of the 
sample may happen, we hope that the good thermal con- 
ductivity of graphene [32j could dissipate the heat fast 
enough. In addition, several techniques such as the use of 
light pulses instead of continuum AC sources might also 
help in the dissipation process. In addition, we have also 
demonstrated that for linearly polarized fields, a change 
in frequency can drive a merging transition (Figjsl), al- 
though it does not necessarily correspond to a transition 
between an I and a SM phase. The resulting new SM 
phase with four Dirac points can clearly be distinguished 
from the one with two Dirac points by a Hall conduc- 
tance measurement, since each valley brings a contribu- 
tion e^ /h{l/2 + n) to the total transverse conductivity, 
where n is the number of Landau levels. 

Finally, for the low frequency regime, where a; <C t, 
the second term of Eq. [T] is very small in comparison 
with the hopping terms and can therefore be neglected. 
According to the symmetry of the couplings in EqJ4] 
for (/:? = 0, a gap at a Dirac point cannot be opened 
by the coupling between the Floquet side-bands [35^. 
Thus, noticing that the renormalization of the hoppings 
does not depend on cj, the position of a PDP can still 
be controlled, and we expect that the PDP will move 
and merge as the field amplitude is modified in low 
frequency regime as well. However, in this limit, a large 
amount of Floquet bands overlap, leading to many band 
crossings for small variations of the amplitude and thus 
to cumbersome quasi-energies behaviors and associated 
phase diagrams. 

Conclusions - In summary, we have shown that the 
merging of the Dirac points in graphene can be driven by 
an AC electric field in a controllable way. We have ob- 
tained rich phase diagrams in which distinct topological 
SM and I phases emerge. At high frequency, uj ^ t^ the 
behavior is easily described by an effective Floquet opera- 
tor, equivalent to the Hamiltonian of undriven graphene, 
but with renormalized hoppings tuned by the field am- 
plitude and polarization. The resulting phases refiect 
interesting topological features, and the critical lines in 
which pairs of Dirac points are created and annihilated 
simultaneously show unexpected localization properties 
in which electrons behave relativistically in one direc- 
tion, while they are localized in the other one. We also 
proved that out of the high frequency regime {uj ~ t), 
and for linearly polarized fields, several PDP can be cre- 
ated, which also involves a merging process which could 
be meassured. At lower frequencies (cj <C t), the descrip- 
tion of the system becomes more involved, but according 
to symmetry arguments we expect that the merging of 



PDP should still be achievable. 

We hope this work could provide an experimentally 
realistic way to control the motion, the merging and the 
amount of Dirac points in graphene. 
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(According to Eq. [s] the crossing of three merging lines 
only is not possible!) 

According to Eq. [5] non linear field polarizations induce 
complex hoppings, and thus break time reversal sym- 
metry. As a consequence, by virtue of the Wigner-Von 



Neumann theorem, the Floquet operator belongs to the shown using an effective Hamiltonian for the side-band 

unitary class, and three parameters (instead of two) are a = when it is coupled to a = ±1. Then, because the 

now required to be varied to find out an accidental de- renormalization does not depend on the frequency, the 

generacy, i.e. a Dirac Point. PDP will move and merge as in high frequency. 
[35] According to EqE] the couplings between Floquet bands 
cannot induce a mass term as far as cj ^ t. This can be 
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Appendix A: S.l Floquet-Bloch theory 

Graphene driven by an in plane AC electric field possesses lattice and time translation invariance: 
H {-x-\- Ri^r -\-T) = JY (x + a^, r) = H {:si^ r -\- T) ^ being a^ the lattice vectors and T = 27r /uj a period of the driving 
field. Then, one can assume solutions in Floquet-Bloch form |^a,k (x, r)) = e'^^'^~'^^'^^^^\ua,\i (x,r)), where e^^k is the 
quasi-energy for the a Floquet state, and k the wave-vector. The Floquet-Bloch states lii^^k (x, r)) are periodic in 
both X and r. The Floquet-Bloch states fulfill the Floquet eigenvalue equation: 

n{k,r)\ua,\c) = ea,k|^a,k), (Al) 

H(k,r) = e-'^-"" {H (r) - idr) e^^-"" 

= H^ (r) - idr (A2) 

where H (k, r) is the Floquet operator. The time dependent Hamiltonian can be obtained through the Peierls substi- 
tution in the tight-binding model, leading to the time dependent hoppings tj (r) = tjC^^^^^'^^ ^ and to the AC driven 
Hamiltonian: 



^k (r) = ^ 5] ^ c {r)l^ c (r),,, tf^ (r) e^^''^ , (A3) 

a, (3 k j 

where we have included the band indexes a, /3, and the time dependence in the annihilation and creation operators 
of Floquet-Bloch states. Expanding in Fourier series the time dependence of c (r)^ j^, and c (r)J^ j^ we obtain: 

^k (r) = E E E E 'T' (^) e"^^(^'"")e'''^4.k.p'C,.k.p (A4) 

a, (3 k p,p' j 

Eq |A4| gives a description of the time dependent Hamiltonian in terms of the time independent operators 
\ Ca,k,p5 <^L k p I • "^^^ quasi-energies are then obtained by direct diagonalization of the matrix: 

((^a,k,p' |H (k, r) |ii/3,k,p)) = ?p/',p - P^Sp'^p6a,(3, (A5) 

rT 3 

(A6) 



?'' P ^ ^ / E e^-<''-^)e^^-^Hf^ (r) dr. 



where we have used the composed scalar product in Sambe space to get rid of the time dependence, and puSp^^p is 
the Fourier space representation of —idt- 



Appendix B: S.2 Ac driven graphene 

The crystallographic structure of graphene consists in a primitive unit cell with two inequivalent atoms (A,B) and 
unit cell translation vectors a2 = f (3,^/3), ai = | (3, — >/3), and as = a (0,0), being a the distance between the 
atom A and B. Each A atom has three nearest neighbors (B type) at positions da = a (1,0), d2 = | (—1, —a/3) and 
di = I (— 1,>/3). In k-space the energy spectrum shows two inequivalent gapless points at K = |^ (l , 4= j and 



K' = 1^ I 1, — ^ 1 called Dirac points. We also define the Time Reversal Invariant Momentum (TRIM) points as: 

Mo = (0,0), (Bl) 

Ms = I (1,0). 
In the absence of driving, the tight-binding Hamiltonian up to nearest neighbors (NN) on the honeycomb lattice is: 

^ (^) = ( p (I) ' V ) . Z' (k) = E he^-"-' > (B2) 

where tj is referred to the hopping along the dj direction. Explicitly, p (k) = X)j=i Pj C^) is given by: 

pi(k) = i, 

P2(k) = ^-'^(3'=^+^'=''), 

P3(k) = ie-''^(3'=--^'=«). 

The spectrum of H (k) has two bands with dispersion relation given by: 



E (k)^ = ±|p (k) I = ±tv/3 + /(k), (B3) 



/ (k) = 2 cos ( avSky ) + 2 cos — /c^^ 



v3a \ /3a v3( 



a , 



/c^ I + 2 cos I — - /Ca^ H -— /c^ 



When the system is coupled to an AC electric field in the dipolar approximation, one can follow S.l to obtain the time 
dependent Hamiltonian and the expression for the Floquet operator in Fourier space. In general, we assume without 
loss of generality that the vector potential for an homogeneous in-plane electric field is given by: 

A (r) = {Ax sin {ujt) , Ay sin {ujt + (/?), 0) . (B4) 

which represents the case of an elliptically polarized field in the x-y plane. Note that cp = reduces Eq |B4| to the case 



linear polarization, and A^ = Ay with (p = 7r/2 reduces Eq B4 to the case of circular polarization. 



To obtain the quasi-energies, one needs to calculate the matrix elements of the Floquet operator by means of Eq |A5| 
and Eq |A6[ The integrals are of the form: 



^ f e^"K^'-^)e^"^^^^'^")e^^^^^^""+^)^r = e^P-pO-^^-Kc+^^cit)) J^,_^ (^a^ + /^^ + 2a/3cos ((^)) . 



r 

T 
which has been solved using the formula: 



e-"'^JAr)= Y. Jv+m(a)J^{P)e- 



m= — oo 



being F = y^^c^I^^^^^^^o^S^osT^, and ^ obtained from the equations: 

{Fcos(^) = a - Pcos{(p) 
Fsin(^)= Psm{(p) 

i.e., 

(3 sin {(f) 



^ = arctan 



a — (3sm{(p) 



Thus, the Fourier space representation of the time dependent tight binding is given by blocks {q = p' — p): 

fC^.P _rr ^^.^_ ( p, (k) 

with Pg(k) = Xlj ^j^g^*^^^ 5 ^1 = a (— 3, >/3) /2, a2 = a (— 3, — >/3) /2, aa = (0,0) and the renormahzed hopping 
elements given by: 



tl^ = tJ,{Az)., 



The functions ^ j and Aj encode all the information regarding the ac field configuration: 



2 



^2,1 = Va2+3A2±2V3A,A^cos((^), 



^2,1 = arctan ^ , ^^ ^^\ ^ , ^3 = ^.a, ^3 = 0. 

\Aj:± v3A^ cos ((/?) y 

In the high frequency regime {uo > 6t), the Floquet bands are decoupled {1-Ln,m is block diagonal), and the system 
can be described by an effective time independent Floquet operator, which is given by the 2x2 matrix: 

^(0) _J e-^^-^ Jo (A) + e-^^-^ Jo (^2) + Jo (^3) A 

V e'^^'^o (^1) + e^^-^ Jo {A2) + Jo (^3) ) ' 



a. High frequency. Different field configurations 
We consider in detail three main configurations for the AC field: 

1. We fix Ax = Ay^/S and vary Ay and cp independently. In this configuration the relation A^ = Ay^/S for the 
AC field amplitudes compensate the asymmetry in x-y of the honeycomb lattice. It gives rise to very simple 
expressions for the renormahzed hoppings. 

2. We fix (/:) = 7r/2 and vary Ax and Ay independently. Note that this is reduced to circular polarization if Ax = Ay. 

3. We fix (/p = 0, and vary Ax and Ay independently. This case contains all varieties of configurations of linearly 
polarized fields. 

The functions Aj and ^^ for the condition Ay = Ax/V^ are given by: 

^2,1 = — ^, ^3 = Axa, 

( sin ((p) \ ^ 
^2,1 = arctan ^^ , ^3 = 0. 

V [1 ± COS (if)] J 

It leads to the renormahzed hoppings t^ = tJo{Axa)^ tf^ = tJo (aA^^ cos (^)), and tf = tJo (aA^; sin (^)). The 
corresponding phase diagram is shown in FigJ6J 
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Figure 6: Phase diagram as a function of the external ac electric field parameters (^ and Ey for the condition Ax — V?>Ay. The 
different colors label the four inequivalent Mj points for the merging. 

The color code in given by: 



Renormalization 


Merging point 


tf + i^+tf = 


Mo (Green) 


^1 ^^ ^2 "•" ^3 


Ml (Brown) 


4 = t( + q 


M2 (Red) 


tr=tr+tf 


M3 (Blue) 



Note that this field configuration allows for the merging of one or two pairs of Dirac points simultaneously (Figj6|. 
For the case of phase difference (f — 7r/2, the functions are given by: 



A 



2,1 



2 v ^ ^ 2/' 



^21= arctan 



A-x 



As = A^a, ^3 = 0. 



being its phase diagram plotted in Figj7[ Note that the case of circular polarization is obtained by tracing a line 



3). 



Ay. Interestingly, the case (p = 7r/2 only presents a trivial insulating phase (description in the main text, Fig. 
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Figure 7: Phase diagram as a function of the external AC electric field parameters Ex and Ey for the condition cp = 7r/2. The 
different colors label the four inequivalent Mj points for the merging. Note that in this case the merging lines Mi and M2 
overlap, giving rise to the existence of two and four simultaneous crossings. 
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For linear polarization [ip = 0) the functions are given by: 



Ai = ^^Al^3Al±2V3A,Ay, 

^2,1 = 0, ^3 = A:,a, ^3 = 0. 



The corresponding phase diagram is plotted in FigjSJ As we can see, this phase diagram has analogies with the one 
in Figj6] 



^ 




Figure 8: Merging lines as a function of the external AC electric field parameters Ex and Ey for the condition (^ = 0. The 
different colors label the four inequivalent Mj points for the merging. 



In addition, in this plot it is clearly seen the imprint of the lattice symmetry (hexagonal pattern) around the 
SMo phase (corresponds to E^ = Ey = of undriven graphene), as we should expect for (f = 0. The effect can be 
understood in terms of the hoppings. By tuning the electric field orientation, we favor a different hopping, leading to 
three different merging lines: blue, red and brown. 



Appendix C: S.3 Lower frequencies 



By decreasing the frequency, the different Floquet bands couple one to each other and the isolated Floquet band 
picture is not accurate anymore (as it was in the high frequency regime, i.e., a; > 6t for graphene). In what follows we 
shall focus in the linearly polarized case ((/? = 0), where the system is time reversal invariant and the mergings can 
happen as in high frequency. Note that, according to Fig|8) its phase diagram contains all possible phases. 

The coupling between different Floquet side bands allows for new types of hoppings, in addition to the spatial 
ones, which include photon absorption and emission processes. These new processes increase the possible transitions 
from the valence to the conduction band, now mediated by photons. This is reflected in the band structure by the 
appearance of a new PDF related with these photon mediated transitions (see FigM. 
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Figure 9: We show the appearance of a PDF at the Mo point when the frequency is decreased from the high frequency hmit. 
We fix Ay = 3.6, Ax — 0, and cp — such that the system is in the NI phase in high frequency. Left plot shows the FEZ at high 
frequency (cj = lOt) in the normal insulating phase. Center plot shows the gap closed by the decrease in frequency (cj = 2.5t), 
where the dispersion relation is linear in one direction and quadratic in the perpendicular one, as it is expected for the merging 
process. Right plot shows the appearance of a new FDF when the frequency is further decreased (cj = 2. It). 



Interestingly we find that a further decrease of the frequency allows one to increase the number of PDF within the 
FEZ (FigflQ|), where the new pairs are associated with photon transitions at the M3 point. 




Figure 10: Creation of a new PDF at M3 in addition to the initial ones by a further decrease of the frequency. Left plot shows 
the case where initially a single PDF is in the FEZ (uj — 1.7). Center plots shows the case uj — 1.6, where the gap is closing at 
M3. Right plots shows the case cj = 1.5, where a PDF has appeared. 



We conclude that in the low frequency regime, and for linearly polarized fields ((/p = 0), the merging is still possible, 

and that a change in the frequency can also be used to produce the merging. For very low frequencies, where cj <C t, 

this mechanisms should remain, although small variations in uj or the field amplitudes would produce significant 

changes in the band structure. Also the existence of an insulating phase is very unlikely for the low frequency, and 

even though the merging would happen, the presence of other PDP would not show insulating behavior (unless the 

field amplitudes are small, such that Aj ex E/uj <C 1 and mainly the zeroth order Eessel function contributes to the 

hopping renormalization, decoupling then the sidebands). 
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